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Abstract 

We investigate the possibility that the late acceleration observed in the rate of expansion of the 
universe is due to vacuum quantum effects arising in curved spacetime. The theoretical basis of the 
vacuum cold dark matter (VCDM), or vacuum metamorphosis, cosmological model of Parker and 
Raval is revisited and improved. We show, by means of a manifestly nonperturbative approach, 
how the infrared behavior of the propagator (related to the large-time asymptotic form of the heat 
kernel) of a free scalar field in curved spacetime leads to nonperturbative terms in the effective 
action similar to those appearing in the earlier version of the VCDM model. The asymptotic 
form that we adopt for the propagator or heat kernel at large proper time s is motivated by, and 
consistent with, particular cases where the heat kernel has been calculated exactly, namely in de 
Sitter spacetime, in the Einstein static universe, and in the linearly-expanding spatially-flat FRW 
universe. This large-s asymptotic form generalizes somewhat the one suggested by the Gaussian 
approximation and the i?-summed form of the propagator that earlier served as a theoretical basis 
for the VCDM model. The vacuum expectation value for the energy-momentum tensor of the free 
scalar field, obtained through variation of the effective action, exhibits a resonance effect when 
the scalar curvature R of the spacetime reaches a particular value related to the mass of the field. 
Modeling our universe by an FRW spacetime filled with classical matter and radiation, we show 
that the back reaction caused by this resonance drives the universe through a transition to an 
accelerating expansion phase, very much in the same way as originally proposed by Parker and 
Raval. Our analysis includes higher derivatives that were neglected in the earlier analysis, and 
takes into account the possible runaway solutions that can follow from these higher-derivative 
terms. We find that the runaway solutions do not occur if the universe was described by the usual 
classical FRW solution prior to the growth of vacuum energy-density and negative pressure (i.e., 
vacuum metamorphosis) that causes the transition to an accelerating expansion of the universe in 
this theory. 

PACS numbers: 04.62.+V, 98.80.Cq, 98.80.Es 
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I. INTRODUCTION 



Recent observations of type la supernovae (SNe-Ia) appear to support the view that the 
expansion of the universe started speeding up about 5 Gyr before the present time, and that 
it was slowing down prior to that time LJ-^. Together with the evidence from the power 

fi n 

spectrum of the cosmic microwave background radiation (CMBR) j5|-jl2| estabhshing that 
the universe is nearly spatially flat, these observations imply that a significant fraction of the 
energy density of the universe is in the form of a dark energy that exerts negative pressure. 
The focus of this paper will be on one possible candidate for this dark energy. 



It has been suggested |13|, that the recent acceleration of the universe may be caused 



by the vacuum energy and pressure of a quantized scalar field in curved spacetime. This 
scalar field is assumed to be free, meaning that it has no non-gravitational interactions. 
Furthermore, the mass m of the field is assumed to have a very small non-zero value, so the 
Compton wavelength of the particle is a significant fraction of the present Hubble radius of 
the universe. The origin of the acceleration in this theory is not the zero-point energy that 
would give rise to a cosmological constant A. The latter would be affected by renormalization 
and would be expected to be very large, if it is not forced to be zero or nearly zero by some 
symmetry, dynamical process, or other principle that is not well understood. The curved 
spacetime effect considered here has a magnitude that depends on the ratio of the mass of 
the free scalar field to a dimensionless constant of order unit y. For a free scalar field in 
curved spacetime, the mass of the field is not renormalized [ifil llfil|. Hence, the mass of this 
free field may be taken to be very small and will remain small within the context of quantum 
field theory in curved spacetime. Thus, the theory studied here, with a very small value for 
m, is internally consistent. 

The gravitational field is treated as a classical field here. One may ask if quantizing 
the gravitational field, by regarding gravitons as quantized fluctuations of the metric on the 
curved spacetime background, would necessarily force the effective mass m of the scalar field 
to be large? The self-energy contribution to the scalar field propagator caused by virtual 
gravitons is evidently not renormalizable, so a cut-off must be introduced. The magnitude 
of the effective mass of the scalar particle would then depend on the cut-off. If the cut-off 
is somewhat smaller than the Planck mass, then the graviton contribution to the effective 
mass would appear to be small. This would not affect our cosmological predictions because 
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they depend on a quantity m^, of order m^, that is determined by the cosmological data. 
We make no attempt here to obtain the value of m from first principles. 

The VCDM model gives a satisfactory fit to the observed CMBR power spectrum and 
SNe-Ia data 17]. It may also lead to the observed suppression of the CMBR power spectrum 
at very low values of / 18[. 

The question of how the very small mass scale fits in with fundamental theories of elemen- 
tary particles and strings is not one that has a clear answer. It has been suggested that 
a scalar particle of the required small mass may arise via symmetry breaking as a pseudo- 
Nambu-Goldstone boson (pngb). It seems possible that such a pngb, when quantized in 
curved spacetime, could lead to similar effects in curved spacetime through an effective 



action j20| or other methods 2l|. In the present paper, only the free field will be considered. 

In the previous work of Parker and Raval Q|, the approximation was made 
that the propagator (the analytic continuation of the heat kernel) was proportional to 
exp [-i (m^ + - 1/6) R) s] (l + {isff2), where m is the mass of the free scalar field, R 
is the scalar curvature, is a dimensionless coupling constant that appears in the equa- 
tion governing the free scalar field in curved spacetime, and /2 is a quantity constructed 
from covariant derivatives of R and contractions of products of two Riemann tensors [see 
Eq. dUD]. 

This approximation was obtained from the expansion of the heat kernel in powers of the 
proper-time parameter s, with the exponential factor involving R coming from summing all 



Alternatively, the 



terms that involve one or more factors of R to all orders in s 
exponential involving R can be obtained by means of a Feynman path integral solution of 
the Schrodinger equation for the propagator 2^ . 

The presence of the exponential in R was shown by Parker and Raval to lead to a growth 
in (T^j/), the expectation value of the energy- momentum tensor of the scalar field, when R 
falls to a value of the order of m?. The reaction back of this growth of {T^^) in the semi- 
classical Einstein equations was shown to cause the expansion of the universe to accelerate 
in such a way as to keep the scalar curvature, i?, of the spacetime nearly constant. 

There are several natural questions that one can raise about the previous approach. For 
example, one may question the validity of the above approximation to the heat kernel, and in 
particular the use of that form of the heat kernel for large values of s. It is the large-s form of 
the propagator that gives rise to the terms in {T^^) that grow large as R approaches a critical 
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value of order m^. A partial justification is that the factor exp [—i {m? + — 1/6) i?) s] in 
the propagator comes from all powers of s, and is thus nonperturbative. However, the series 
in powers of s is in general asymptotic, not convergent, so the convergent subset of terms 
involving R that are summed may not fully reflect the form of the propagator at large s. It 
should be mentioned that the Gaussian approximation to the path integral of Bekenstein and 



Parker 



2J| is independent of the power series in s and gives reason to expect an exponential 



factor at large s. However, it is also an approximation, and thus leaves room for some 
modification. 

In the present paper, we directly examine the large-s asymptotic form of the heat kernel 
for several cases in which the heat kernel is exactly known. We postulate a general expression 
for the asymptotic form that is consistent with the cases studied. We show that the transition 
to constant R occurs as a result of the large-s asymptotic form of the heat kernel. In our 
numerical integration of the Einstein equations in the FRW universe, we use the postulated 
asymptotic form with a particular choice of a factor that is quadratic in the Riemann tensor. 
This choice (namely, f2) is consistent with the exact asymptotic forms considered, but 
there are also other consistent possibilities, as we discuss. We expect that the transition 
is fairly generic, as the choice we make of the quadratic factor is sufficiently complicated 
to be representative. As noted below, the possible runaway solutions of the field equations 
obtained from this higher derivative effective action do not occur with physically acceptable 
initial conditions. This approach based on the asymptotic form of the heat kernel is general 
enough that it may be useful in other fundamental theories, such as string theory, that have 
higher derivative terms in their low-energy effective actions. 

One may also question the validity, in the earlier work of Parker and Raval, of neglecting 
derivatives of the Riemann tensor in arriving at the expression for (T^j^). The reason for 
neglecting such derivatives was that the present expansion of the universe is slow enough 
that terms involving derivatives of invariants were expected to be small with respect to other 
terms in the effective action that do not involve such derivatives. However, the solution for 
the expansion of the universe that they found went through a relatively rapid transition from 
a matter-dominated expansion to a constant-i? expansion. During that transition it may be 
necessary to include the derivative terms that were neglected. Furthermore, by including 
the higher derivative terms one may generate unstable runaway solutions of the Einstein 
equations. 
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In the present work, we keep the derivative terms that were neglected in the earher 
work. Then the covariant conservation of {T^^,) is satisfied at all times without neglecting 
derivatives. We find that the transition to constant R still occurs. We also carefully consider 
runaway solutions. We find that if the early evolution, prior to the transition, is described 
by the usual classical solution to the Einstein equations (as one expects because {T^J) is 
negligible before the transition), then there is a transition to the constant R solution with 
no runaway solutions. There are runaway solutions only if one takes the initial evolution to 
be significantly different from that of the classical solution prior to the transition. 

We also generalize the previous work by including the neglected effect of dissipative 
processes. These could come from very small interactions of the scalar field with other 
fields, which could involve dissipative processes such as very slow decay of the scalar particles 
into less massive fields, or weak radiative processes, such as the production of gravitational 
waves. In the present paper, we model such dissipation phenomenologically by introducing 
a small, but nonzero, imaginary part in the mass term. 

In their previous work, Parker and Raval mentioned the possibility that the mechanism 
they considered may be relevant to early infiation. However, because the series in powers 
of s is more difficult to justify as an approximation when the curvature terms are large, 
they applied their work only to the recent universe. However, now that the mechanism 
has been further justified by considering known exact asymptotic forms of the heat kernel, 
the time may be ripe to apply this mechanism to large-mass scalar fields that are present in 
elementary particle theories. The acceleration effect will still come from the large-s behavior 
of the heat kernel using the same postulated asymptotic form. The fact that the mass is 
large would not seem to prevent the growth of {T^y). Because of the large mass, these 
bosons could give rise to early infiation. At the very least, our mechanism could modify the 
considerations of early infiation that come from possible self-interaction potentials of these 
fields. It is also possible that our mechanism could give rise to a satisfactory new model of 
early infiation in the absence of self-interaction potentials. Thus, this type of field might 
serve as a candidate for an inflaton. The dissipative term would then be larger because of 
the interactions and decay channels of the massive scalar boson. One would of course have 
to consider the questions of reheating and the perturbation spectrum produced in such an 
infiationary epoch. This mechanism for early infiation is worthy of investigation, but it is 
not our main focus here, and will not be considered further in the present paper. 
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The very low-mass, free, scalar particle that we consider in this paper as a possible cause 
of the recent acceleration of the universe, may appropriately be called an acceletron, to dis- 
tinguish it from an inflaton. Even in the absence of other mechanisms for early inflation, the 
acceletron^tself would produce an early inflation at the Planck scale, similar to Starobinsky 
inflation |25^. Assuming that a successful exit from early inflation occurs, then at a much 
later time, the same acceletron would produce the presently observed acceleration of the 
universe. 

In Section |TT1 we outline the basic theoretical structure, including the regularized ex- 
pression for the effective action. In Section IIIH we discuss the asymptotic form of the heat 
kernel for large-s, and make a postulate for this asymptotic form that is consistent with the 
asymptotic forms of several exact solutions for the heat kernel. In Section HVl we show that 
the large-s asymptotic form of the heat kernel is what gives rise to the terms that make large 
contributions to the vacuum energy- density and corresponding negative pressure when the 
scalar curvature falls to a magnitude comparable to the square of the mass of the acceletron 
field, (f). In Section |3 we write the renormalized effective action, and give the corresponding 
energy-momentum tensor. The relevant variations, including terms with derivatives of the 
Riemann tensor, are given in Appendix ^ In Section IVIl we briefly summarize results of 
numerical integration of the Einstein equations in an FRW cosmological spacetime. More 
detailed results of the numerical integration will be given in a later paper that is now in 
preparation. In Section FVIH we present our conclusions. 

II. EFFECTIVE- ACTION FORMALISM 

In this section we outline the method of relating the one-loop effective action of a quantum 
field in curved spacetime to the generalized ^-function and the propagator in the proper-time 
formalism. 

The action S of an uncharged scalar field (j) in curved spacetime can be written in the 
form 




(2.1) 



with H{x) being the differential operator 




(2.2) 
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□ = g^^''V and ^ being a coupling constant between the field and the scalar curvature 
R = R{x) of the spacetime. 

The corrections to the classical action of the gravitational field that are caused by quan- 
tum fluctuations of (p can be found by evaluating the effective action Wg = Wq[g^^]. The 
latter is defined for a given spacetime with metric g^i, through the functional integral 

e^^« = j d[(j)] e'^. (2.3) 

In analogy with the Feynman path integral, the functional integral is proportional to a prob- 
ability amplitude to go from an "in" state to an "out" state determined by the configurations 
of (f) that one sums over. Rather than explicitly specifying the initial and final configurations 
of (f), the "in" and "out" states are specified implicitly by boundary conditions used later 
in evaluating the functional integral. [To insure convergence of the functional integral, the 
operator H{x) is understood to be H{x) — ie, where e is a positive infinitesimal.] 

A ,.o„^ro i,„agi„a,v part of W, implies a — rate of part.de producfo,. of this 
scalar field |2a|. If the spacetime is such that we can neglect the imaginary part of Wq, then 
the "in" and "out" states of the field are equivalent, differing only by a phase factor. The 
expectation value of the energy-momentum tensor of the scalar field in this state is then 
given by 

= -L^. (2.4) 

In the semi-classical theory, this expectation value appears on the right-hand-side of the 
Einstein gravitational field equations, along with the other sources of the gravitational field. 
Later, we will apply this formulation to a spacetime for which it appears quite suitable; 
namely, the Friedmann-Robertson- Walker (FRW) universe describing the averaged behavior 
of the recent expansion of the universe. 

Because the action S of the free scalar field is quadratic in 0, the functional integral 
reduces to a Gaussian integral which gives (to within a constant normalization factor that 
does not affect the variation of Wq): 



l/^Det(^//.2), (2.5) 

where H is the abstract operator defined through 

{x\H \x') = H{x)5{x,x') , (2.6) 



with \x) being position eigenstates normalized through 



{x\x') = 5{x, x') = ^'(^^') . (2.7) 



Then 



Wq = (V2) InDet (////i^) = (i/2)Trln (^H/jJ^^ , (2.8) 

with fi an arbitrary constant having the units of mass. The renormahzed effective action 
will not depend on /i. 

The expression for Wq can be regularized by writing it in terms of the generalized (- 



function 



27 



28|, 



We can rewrite Eq. ()2.8|) as 



. The latter is defined as 

C(z/) = Tri^'-'^ = Tre"^'"^-^. (2.9) 



Wq = ~l[aO) + Hf^')aO)], (2.10) 



where ('{i^) = dC,{v)/dv. This expression can be regularized by analytic continuation of C,{v) 
in the parameter v to make it and its first derivative well-defined at = 0. The dependence 
on the arbitrary constant /i can then be absorbed into the definition of the renormalized 
constants such as Gat and A that appear in the Einstein action. As is well known, additional 
terms quadratic in the Riemann tensor must be added to the Einstein action to absorb all 
the dependence on /i through renormalization. 

It is convenient to introduce an integral representation of the generalized (^-function: 

C(z/) = Tr^-- 

= 1i\v{u)-^ idsiisy-^-''^''"''^^. (2.11) 

The operator H{x) of Eq. ()2.2|) may be regarded as the Hamiltonian of a fictitious nonrela- 
tivistic particle moving on a curved 4-dimensional hypersurface having coordinates x^. The 
operator exp(— isif) is the quantum mechanical evolution operator of this fictitious particle, 
with s being the "proper time." This idea was introduced by Schwinger in the context of 
quantum electrodynamics and applied by DeWitt in curved spacetime. The trace can be 
represented as an integral over a complete set of position eigenstates |x) of the fictitious 
particle: 

(^(u) = [ d^xy/^ [ ids{isy-\x\e-'-<^-''^\x). (2.12) 
^{^) J Jo 



The propagation amplitude for the particle to go from position x' at time s = to position 
X at time s is 

K[x,x']is) = (x \e-^<^-^')\ x') . (2.13) 

This propagator satisfies the Schrodinger equation 

d 

i^rKix, x'] is) = H(x)K(x, x'; is), (2.14) 
as 

and satisfies the "initial" condition 

\imK{x,x';is) = 6{x,x'). (2.15) 

Note that what appears in Eq. ()2.12|1 is the coincidence limit x' ^ x oi Eq. ()2.13j) . If one 
were to replace is by s in the Schrodinger equation, then it would become the equation 
governing heat flow on the 4-dimensional hypersurface, in which case K{x, x', s) is the heat 
kernel. Thus, the heat kernel and propagator are related by analytic continuation in the 
proper time s. 

The asymptotic expansion of the propagator in powers of s can be obtained from the 
Schrodinger equation by iteration. The proper time s has the dimensions of length squared, 
so the coefficients of successive powers of s in the series contain contracted products of 
increasing numbers of curvature tensors and/or covariant derivatives to balance the dimen- 
sions. The leading terms in this power series covariantly characterize the short wavelength 
behavior of the quantum field. In four dimensions, the terms of order s~^, s~^, and s° in 
the series are the ones that would give rise to ultraviolet (UV) divergences in the unregu- 
larized expression for the expectation value of the energy-momentum tensor. Those are the 
terms that are absorbed through renormalization of the coupling constants of the curvature 
terms in the classical Einstein action (with counterterms quadratic in the Riemann curvature 
tensor included). 

Writing the coincidence limit, K{x,x]is), of the propagator (in four dimensions) in the 
form 

K{x, X- is) = -r4r^ e-^(™'--)^F(x, x; ts), (2.16) 
the proper-time series for F{x, x; is) has the form 

oo 

F{x,x;ts) = Y,{^s)^fj, (2.17) 

j=0 

10 



where the first three terms are given by 



/o = fo{x,x) = 1, 
/i = fi{x,x) = -^R 
f2 = f2{x,x) = ^eR 



(2.18) 
(2.19) 



^ ' ^ [(1 - 30e) OR + R^p^sR"^^^ - RapR""] , (2.20) 



2 ' 180 

with ^ = ^ — 1/6. In general, this proper time series for F{x,x;is) is an asymptotic series 
that does not converge to the exact solution for F{x, x; is). Thus, it can only supply limited 
information about the large s behavior of F. 

It turns out that the Schwinger-deWitt proper-time series given in Eq. ()2.17|1 is sufficient 
to calculate one of the terms appearing in Eq. ()2.10|) . namely the term involving C(0)- 
order to do so, we first have to analytic continue Eq. ()2.12|) to the value u = 0. Substituting 
Eq. ()2.16|1 into Eq. ()2.12|) and performing three integration by parts, one has: 

i 



a^) = - 



167r2(z/-2)(z/-l)r(i^ + l 

^3 



d Xy/—g 



'^"''-''^'F{x,x;is) 



(2.21) 



which is regular at z/ = 0. In fact, one easily obtains 



C(o) 



327r2 
i 

32^ 



d Xy 
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d{isy 

j d'xy/^g K - 2mVi + 2/2) . 



s=0 



(2.22) 



The remaining term in Eq. (j2.10p involving C'(0)) the other hand, cannot be calculated 
only based on the small-s behavior of F{x,x;is). Taking the derivative of Eq. (j2.2H) with 
respect to u, and applying the result to z/ = 0, we have 



C'(0) 



5 + -')«<»>-3^ 



X 



ids In(zs) 



93 



d'^Xy/^ 



<rn^-i-)^F{x,X;is) 



(2.23) 



'0 d{is)' 

where 7 = — (d/dz) lnr{z)\^^-^ is the Euler's constant. As we shall show in Sec. IIVI a most 
important contribution to C'(0) comes from the large-s regime of the integrand of Eq. ()2.23p . 
Combining Eqs. ()2.22j) and ()2.23|) into Eq. ()2.10|) . we have the following expression for the 
one-loop effective action: 



~6^I ^'^v^ t^^O' - i"^' - + 2/2) In ft'] , (2.24) 
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where, for given < a < (3 < +00 



/3 ^3 

ids ln(2s)— — 

o[is) 



^-i{ni^-ie)s p/ 



(2.25) 



and jl is such that In/i^ = 3/2 + 7 + In /i^. 

Before analyzing in detail the properties of X(0, +00), we turn to some known results 
about the coincidence limit of the propagator, K{x^x;is). In the next section, we consider 
several cases in which the exact solution for K{x,x\is) is known. Based on these cases, we 
postulate a general form for the large-s asymptotic form of K{x, x; is). A similar asymptotic 
form is also suggested by the summation of the subset of terms in the proper-time series for 



22 , . This summation 



F{x,x; is) that involve one or more factors of the scalar curvature R 
involves all powers of s, and thus may give some information about the asymptotic form of 
F{x,x;is) for large s. In addition, the Feynman path integral solution of the Schrodinger 
equation ()2.14|) to Gaussian order 2J] suggests a similar large-s asymptotic behavior. It 
is the asymptotic form of F{x,x;is) for large s that determines the contributions of long 
wavelength infrared (IR) quantum fluctuations of the field to the functional integral for the 
effective action. If the mass of the scalar particle is very small, it has a very large Compton 
wavelength which may significantly influence the long wavelength quantum fluctuations of 
the (f) field. As we shall see, there are spacetimes in which this influence can cause the 
expectation value of the energy-momentum tensor of the free scalar field to become large. 



III. IN SEARCH OF THE LARGE-5 BEHAVIOR OF THE PROPAGATOR 

It is well known that the series for F{x, x; is) given in Eq. ()2.17|) is not convergent 



in general. However, as conjectured by Parker and Toms [22|, and later proved by Jack 
and Parker Q], Eq. ()2.17|1 contains a convergent sub-series involving powers of the scalar 
curvature R, with the property that when it is summed, the (asymptotic) series that is left 
does not possess any power of the scalar curvature (without derivatives applied to it). More 
precisely, the series for F{x, x; is) in powers of s can be written (with ^ = ,^ — 1/6) as 



00 



F{x,x;is)=e-'^^^J2^ts)^fv (3.1) 

j=0 
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with the first three terms being 

/o = fo{x,x) = l, (3.2) 
/i = /i(x,x) = 0, (3.3) 
h = Mx, x) = -^ [{1 - 30e) DR + Rap.sR'^''^' - RapR""] , (3.4) 

and the subsequent terms being free of any terms containing undifferentiated factors of R. 
Obviously, fixing a particular background might allow one to write the contractions of the 
Riemann and Ricci tensors appearing in the fj in terms of the scalar curvature R. For 
example, the symmetries of de Sitter spacetime (a case that will be analyzed later) make 
it possible to express any geometrical scalar quantity in terms of R alone. However, the 
dependence of fj on R would be different for different backgrounds. Similarly, there are 
identities that depend on the dimension of spacetime that would allow one to reexpress 
some of the coefficients fj in such a way that terms involving R would appear 
The expressions for the /-,• that are valid for general metrics and in spacetimes of arbitrary 
dimension do not include factors of R. The form of the coincidence limit of the propagator 
given in Eq. ()3.Hl is called the i?-summed, or the partially-summed, form of the propagator 
or heat kernel. The factor of exp{—iC,Rs) in Eq. ()3.H1 sums a covariant and dimensionally- 
invariant set of terms to all orders in s. Therefore, it may contain information about the 
large-s behavior of K{x,x;is). Earlier, Bekenstein and Parker [2^, using Fermi coordinates 
in curved spacetime, were able to obtain the Gaussian approximation to the Feynman path 
integral solution of the Schrodinger equation ()2.14|) . This approximation did not restrict s 
to small values. The result they obtained for K{x, x'; is) reduces in the coincidence limit to 
the result obtained from the first term of the series in Eq. 1)3.111 . This gives further support 
to the impression that the exponential factor has relevance to the nonperturbative large-s 
behavior of K{x,x;is). 



In the work of Parker and Raval |l3l.ll4l| they considered the effects of this nonperturbative 
term in K{x, x; is) on the effective action. Here we clarify and extend their work in a 
number of respects. First, by considering several exactly known heat kernels we generalize 
the expression for the large-s behavior of K{x, x; is) that was suggested by the i?-summed 
form and by the Gaussian approximation. Second, we prove that the growth in the vacuum 
expectation value of T^i, that occurs in certain spacetimes comes directly from the large- 
s asymptotic form of the propagator. Third, we incorporate the possibility of dissipative 
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processes into our expression for the effective action. We also summarize numerical solutions 
of the resulting semi-classical Einstein equations that we obtained in the FRW universe. 
These numerical solutions include the contributions of higher derivatives of the Riemann 
tensor. 

In order to arrive at a sufficiently general conjecture for the large-s asymptotic form of 
K{x,x;is), let us analyze three particular cases where the propagator K{x,x';is) is known 
exactly, namely, de Sitter spacetime, the Einstein static universe, and the linearly-expanding 
spatially-flat FRW universe (the latter one assuming conformal coupling C, = 1/6). The exact 
Euclidean heat kernel Ke{x, x'; s) (which is the Euclideanized form of the propagator K) in 
de Sitter spacetime for a massive scalar field was calculated by Dowker and Critchley 35| 
and given by: 

= w + \i - + I ^^^^^ ' (^-^^ 

where p = cos{\^2a/a), a = a{x,x') is half the square of the geodesic distance between 
X and constant related to the scalar curvature by R = 12/a^, and Pj are the 

Legendre polynomials. [Here we use signature (—,+,+,+), so that we have performed 
the substitution a —a in the results of Ref. |3,9], which uses signature (+,—,—,—).] 
Using dPj{p)/dp\p=i = j{j + l)/2, factoring out the j-independent exponential factor, and 
multiplying Eq. ()3.5p by i to go from the Euclidean to the Lorentzian metric, we obtain the 
exact de Sitter propagator in the coincidence limit: 



-e 



R2 



87r2 144 

°° l-iRs 
x2^j(j + l)(j + l/2)exp 

i=o 



12 



J(J + 1) 



(3.6) 



where = + ^R. (The exponential factor multiplying iKE{l,s) sets the propagator 
to the notation used here.) Notice that the exponential for the term j = would be 
exactly the same as the exponential obtained from the i?-summed form of the propagator, 
g-iAf^s j^Q-^gygr^ ^j-^ig exponential in Eq. ()3.6|) appears multiplying a factor proportional to 
j, and therefore gives no contribution to the summation (possibly due to the high degree 
of symmetry of de Sitter spacetime). Since the effective action Wq depends linearly on 
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the propagator K [see Eqs. ()2.im) . ()2.12p . and ()2.13p ]. we can investigate separately the 
contributions to Wq coming from each term in the summation of Eq. ()3.6p : 

K,{x,x-is) ^ -^^e-^K)-)^7^2a)(^s)^ 3 > 0, (3.7) 

with Mj) = M^+]{] + l)i?/12 and 7^20•) = j{j + l)(j + l/2)i?V72. Note that 

00 

K{x, X ;is) = ''^^ Kj{x, X ;is). (3.8) 
i=o 

As will become clear in the next section, the dominant term in Eq. (3.8) for large s is the 
one with j > for which is the smallest. This value of j, which we denote by k, will of 

course depend on the value of ^. For example, if ^ > —1/6, then k = 1 and M^^^ = + ^R- 
More generally, we can say that (in the distributional sense) the dominant behavior of the 
propagator in de Sitter spacetime, for large s, is given by 

K{x,x-ts) ~ -_L_e-^('^«~^^)^7^2(.)(^.)^ (3.9) 

with k>l. 

Our second example, is the Einstein static universe. In this case, the exact form of the 
coincidence limit of the propagator is given by 



K{x, X, IS) = - ,-^iM^-^^)s <j 1 + 2 5^ e^^-^-^'^V^ (1 + %f^a'/s) \ , (3.10) 

3=1 



where a is related to the scalar curvature through R = Q/a^. The summation in j appearing 
in Eq. ()3.1()|1 is related to the fact that in the Einstein static universe there are infinitely 
many geodesies connecting any point x to itself. The contribution to the propagator given 
by the direct path connecting x to itself (i.e., the trivial path), is encompassed by the factor 
1 inside the curly brackets of Eq. ()3.10|) . Note then that the Gaussian approximation for 
the propagator gives the exact direct-path contribution in the Einstein static universe. If we 
restrict our attention to the direct-path contribution, we have the large-s regime of Eq. ()3.1()|1 
given by 

We also find that approximation of the infinite sum in Eq. ()3.10|) by an integral strongly 
suggests that the contributions of the indirect paths sum to a quantity that grows much 
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more slowly than s^. Therefore, we conclude that the rate of growth of the de Sitter heat 
kernel for large s is faster than that of the heat kernel in the Einstein static universe. 

As our last example, let us consider the spatially-flat FRW universe that is expanding 
linearly in the FRW proper time coordinate. With conformal coupling, ^ = 1/6, between 
the field and the scalar curvature, the propagator was calculated by Chitre and Hartle 



and Charach and Parker 



3^. In the coincidence limit, it can be put into the form 

oo g-6i(coshD)2/(i?s) 




K{x, X- is) = - -—-^ e-'^^ Ji + p dv — . } , (3.12) 



where (3 is an arbitrary constant. Note again that the overall exponential factor appearing 
in Eq. ()3.12|1 is consistent with the i?-summed form of the propagator, recalling that in the 
case of conformal coupling one has ^ = 0, i.e., = m?. Another point worth mentioning 
about Eq. ()3.12|) is that for large values of s, the integral inside curly brackets gives no 
contribution to order for j > 0. [This conclusion can be drawn by taking the limit s oo 
in the integrand in Eq. ()3.12j) and noting that the integral vanishes in this limit.] Therefore, 
we have that in the large-s regime the dominant behavior of Eq. ()3.12|) can be written as 



K{x,x-is) ~ --— ^e-*("^'-*^)^7^A(^s)^ (3.13) 

where A is some number smaller than 1 and TZ\ is a scalar quantity with the same dimension 
as R'^. For ,^ 7^ 1/6, we would expect m? in ()3.13|1 to be replaced by M^. 

Of the FRW metrics considered so far, the large s asymptotic form of the heat kernel of 
de Sitter spacetime, Eq. ()3.9|) . has the highest power of s multiplying the exponentials that 
appear in all the examples. Therefore, in arriving at an ansatz for the large s asymptotic 
form of the heat kernel in a general (but not pathological) FRW universe, we must include an 
exponential multiplied by a power of s that is at least as large as the power that appears in 
()3.9|) . Therefore, based on the results and discussion presented so far in this section, about 
the form of the coincidence limit of the propagator, K{x,x]is), in de Sitter spacetime, the 
Einstein static universe, and the linearly-expanding spatially-fiat FRW universe, it seems 
reasonable to make the ansatz that, at least in the four-dimensional FRW universes that 
will be considered in a later section, the general form of the dominant term in the large-s 
behavior of K{x, x; is) is given by Eq. ()2.16|) with 

F{x, X] is) ~ Un e-'''-^'{is)'' (3.14) 
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for some integer n (to be discussed in the next paragraph), with Xn being a dimensionless 
number and TZn = T^nix) a scalar quantity constructed from the metric g^u and having the 
same dimension as the n-th power of the scalar curvature R. It seems clear that the value 
of Xn will depend on the constant ^ that appears in the field equation for (p. It may also be 
related to the topological properties of the spacetime because the large-s asymptotic form 
of F{x,x;is) could conceivably sense the large-scale structure of the spacetime. We will 
assume that the value of ^ is chosen such that Xn is negative, as negative values of Xn are of 
interest for the cosmological effect that we consider in a later section. The parameter that 
determines the cosmological effect is the ratio —m?/xn = fn^- 

The particular cases analyzed here may give some hints on the values of n and 7?.„. Since 
we are looking for the general dominant large-s behavior of K{x,x;is), the de Sitter case 
analyzed above seems to suggest that n is 2. (It could be larger than 2, but our examples give 
no evidence of that.) Then the Einstein static universe and the linearly expanding universe 
would be special cases in which subdominant asymptotic terms become dominant as a result 
of the vanishing of the coefficient of the dominant asymptotic term. Should this be true, 
then TZ2 must be a scalar that vanishes when calculated in the Einstein static universe and 
in the linearly-expanding spatially-flat FRW universe, while being non-zero when calculated 
in de Sitter spacetime (recalling also that its dimension is the same as the dimension of R^). 
It is not difficult to verify that some of the candidates for the general expression of 7^2 in 
four dimensions are given by the integrand of the Gauss-Bonnet invariant, 

G = Raf^^sR""^^^ - ^RapR""^ + R\ (3. 15) 

the second-order term in the /2-summed form of the Schwinger-deWitt proper-time series 
[see Eq. (ED)], 

f2 = z^ [(1 - SOOOR + R^p.sR'^^^' - R^pR''^] , (3.16) 

ioU 

and the scalar quantity 

S = R^pR^^-^. (3.17) 
Actually, linear combinations of G, /2, and S are also candidates for the general expression 

of n2. 

It is already clear from the work of Parker and Raval, that it is the exponential that 
appears in the asymptotic form of the propagator that is the key to causing an acceleration 
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of the expansion of the universe. In our numerical work, to be discussed in a later section, 
we will take n = 2 and a particular choice for 7^2- We will also mention numerical results 
we obtained by taking n = 0, which confirm that the exponential is responsible for bringing 
about acceleration of the expansion. Although the case n = was worth considering because 
it involves fewer time derivatives of the metric, the case n = 2 appears to behave better, as 
will be discussed. 

It is important to stress that the results we shall derive in the next section are independent 
of the specific form of the term TZn in Eq. ()3.14|) and of the assumption that n = 2 (as long 
as n > and TZn is not identically zero for a general metric, which one immediately sees to 
be true from the cases we have considered). For this reason, in the next section we will use 
the generic ansatz for the large-s behavior of F{x, x; is) given in Eq. (j3.14|) and only in later 
sections will we take n = 2 and assume a particular form for 7^2- 

IV. LARGE-S' ASYMPTOTIC BEHAVIOR OF THE PROPAGATOR AND NON- 
PERTURBATIVE INFRARED QUANTUM EFFECTS 

Returning to the form of the effective action given in Eq. ()2.24|) . we will split the quantity 
X(0, +oo) in two terms [see Eq. ()2.25|) ]: 

J(0,+oo) = J,eg+ J/i?, (4.1) 

where 

^reg = ^T^eg ( A/r) = X(0, A/^j) , (4.2) 

IlR = IlRiXlR)=IiXlR,+(x), (4.3) 

with XiR some "large" but fixed parameter with dimension (length)^. The contribution 
of long wavelength (IR) fluctuations of the field to the effective action is given by the 
integration over large s. 

We will be most interested in the quantity Ijr and will assume is a well-behaved 

function of the metric if x is a non-singular point of the spacetime. This assumption seems 
quite reasonable because we do not expect the function F(x, x; is) appearing in the definition 
of Ireg to givc any problem in the limited interval of integration (0, A/r). Moreover, from 
the Schwinger-deWitt proper-time series we even know that F{x,x;is) — 1 as s — > 0. The 
situation is different, however, for the large-s contribution Ijr, as we will analyze next. 
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Combining Eqs. ()4.3|) . ()2.25p . and ()3.14p . we have, for sufficiently large Xir, 



TiR / ids In (is 



^3 



d{is) 



-i{Ml-ie)s /^\n 



2\n 



{Ml-ief J{\iR,Ml 



where 



(4.4) 



Ml = m' + XnR 



(4.5) 



and the quantity J is defined by 



ids In(is) e 



-i{Ml-ie)s 



~i\lR{Ml~ie)s 



= i^iR J ds[\ias + \ia{iXiR)]e 
ln{iXiR)e-'^''^^^^'-''^ 

/oo 
ds {In 3)6-'^""^^^"-''^% (4.6) 

with s = s/XjR. Using Eq. (4.358.1) of Ref. js^, the integral appearing in Eq. ()4.6|) can be 
evaluated: 



: ^ ^ " ^p\[^X,R{M-^-^eW 



/3=1 



[tXiR{M^ - {7 + ln[a,fi(M2 - te)]} 



j=0 



(4.7) 



where T{/3, a) is the incomplete gamma function and in passing from the first to the second 
line of Eq. ()4.7|) we have used the expansion of r(/5, a) in powers of a (see Eq. (8.354.2) of 
Ref. j3j]). Then, Eq. (gS)) becomes 



J{XiR,M^) = -{M^-^e)-^[^ + \n{M'^~^e 



'-IVHXirM^V 



j=0 



U + 1) 



ln(iA 



IR) 



(4.8) 



Finally, using Eq. ()4.8|) to evaluate Eq. ()4.4p . we obtain for the dominant infrared contribu- 
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tion to the effective action 



QT. 



+(-l)"7^ 



d{Ml) 

& 



d{MlY 



j=3 



(j-2)! 



J 



IR) 



(4.9) 



Notice that the summation now starts at j = 3. 

There are a few points worth mentioning about Eq. ()4.9|1 . First, note that the series 
appearing above is (absolutely) convergent for any (finite) value of (A/ijM^). Moreover, 
the result of any number of derivatives with respect to applied to this series is still 
(absolutely) convergent. These facts imply that the first term on the right-hand-side of 
Eq. ()4.9|) is the dominant one when is sufficiently small. (In particular, for n > 2, the 
value of this dominant term is unbounded when — > 0, i.e., when R —>■ —irP'/xn)- Notice 
that this dominant term in Eq. ()4.9j) is independent of Xir. Thus, the effective action given 
in Eq. ()2.24|1 can be approximated by 

(_l\n+l 



d Xy 



gn 



d{Ml 



2\n 



le 



(4.10) 



where 



q 



^ j d'^Xy/^!^ Ireg - {m^ - 2m^ fi + 2/2) lu /i' 

(7 + lnm=')M^ 



647r 
+ (-l)"7^ 



d" 



+ (a 



IR) 



\-2 



00 

E 

J=3 



-lyitXjRM^y 



\n{iX 



IR, 



(4.11) 



(j-2)! Vj-2 

denotes the part of the effective action whose value and variation with respect to the metric 
Qp^p are well-behaved for any value of M^. 

The positive infinitesimal e in Eq. (j4.10|) was originally introduced to ensure convergence 
of the functional integral for exp(iiyg). Dissipation through small interactions with external 
systems can be modeled by allowing e to be a finite quantity, small with respect to m?. For 
example, a long but finite lifetime of the particle associated with the field may be modeled 
in this way. 
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V. RENORMALIZATION AND VARIATION OF THE EFFECTIVE ACTION 



Now, let us carry out the renormalization of the effective action given in Eq. ()4.1Up . 
Notice that because we have used the (^-function formahsm, the effective action Wq aheady 
at this stage is free from "unphysical" divergences. However, Wg given in Eq. (|4.1U|) depends 
on the unobservable parameter /i through terms up to second order in the curvature of the 
spacetime [see Eq. ()4.1H) ]. Adding a 6are gravitational action containing terms up to second 
order in curvature, 



Wg = j rf^xv^ (-2kA + kR + aiR^ + a2R^uR^'' + a^Rf^^p^R^^'f"') 



(5.1) 



with K, A, «!, 0:2, 03 being bare gravitational constants, we can absorb jl into the definition 
of observable low-curvature gravitational constants Kq, Aq, aio, «2o, «3o with the result that 
the low-curvature limit of the total action has the form 

W = Wg + Wi^ J d^x^^ {-2koAo + KoR + aioR^ 

+a2oR^.uR^'' + asoRf.upaR^""' + ...) . (5.2) 

Since this procedure was performed by Parker and Raval in Ref. [3], using the Schwinger- 
deWitt proper-time series, we will skip this detailed calculation here. What we will compute 
now is the form of the total effective action after renormalization. In order to do so we note 
that the renormalization procedure described above has the net effect of replacing the bare 
gravitational constants in the bare gravitational action by the observable ones (thus giving 
the renormalized gravitational action), while adding to Wq terms up to second order in the 
curvature in such a way to cancel the terms up to second order originally in the low-curvature 
expansion of Wq (in this way giving the renormalized effective action). Therefore, we have: 

w = Wg + Wg = (iy,)_ + (iy,)_ , (5.3) 

with 

(^s)ren = j {-2KoK + f^oR + ai^R^ 

+a2oRp.R^'' + asoR^.^paR'^'n (5.4) 
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and assuming, for example, n = 2 in Eq. ()4.1Up motivated by the discussion of Sec. |^ 

m... - (^r).„-^/ c^^.v-^7^.ln(^) 



rf^xv/^7^2ln(^^^^'). (5.5) 



647r^ J \ 
Recall that M| = + X2-R [see Eq. ()4.5p ] and that 7^.2 is quadratic in the curvature (see 
Sec. In passing from the first to the second line of Eq. (j5.5j) we have used ln(M;^ 
In |M| + it\ + i7re(-M|), where 



2 -«ej ~ 



^ , X < 



1/2 , X = 0, (5.6) 
1 , X > 



is the Heaviside step function. The quantity {^q^^)^^^ in Eq. ()5.5|) is obtained by applying 
the renormalization procedure to WT*^^ alone. 

n 

The imaginary part of Eq. (j5.5p is related to particle production [IJ]. However, since the 
order of magnitude of the imaginary term written explicitly in Eq. (j5.5|) could be comparable 
to the order of magnitude of (W^g^^)^g„; which may also have an imaginary part that we are 
not obtaining explicitly, we will not analyze this phenomenon here. 

Setting the observable gravitational constants aio, a2o, and a^o to zero (to reproduce 
the classical vacuum Einstein equations, with a cosmological constant, in the low-curvature 
limit), the condition that the variation of the total effective action vanishes for arbitrary 
variations of the metric gives us the semi-classical vacuum Einstein equations: 

Rf,u - ^ gi^uR + Kg^lu = SuGn {Tf,^) , (5.7) 

where Gn = l/(167rKo) is Newton's constant, and the vacuum expectation value of the 
energy-momentum tensor of the scalar field is defined through 

As noted earlier, we will set Aq = to see if the expression for {T'^") could explain the 
observed acceleration of the expansion of the universe. 

In order to numerically integrate the semi-classical Einstein equations in an FRW uni- 
verse, we must adopt an explicit form of (T^*^). For this purpose, we assume again n = 2 
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and adopt an expression for 7^2 that satisfies all the properties discussed in Sec. lIIII Namely, 
we let 

n2 = ah = ^ [(1 - m)^R + RaMsR'"^^' - Ro.pR'^] , (5.9) 

where a is some dimensionless constant. This is the term that appeared in the i?-summed 
form of the heat kernel. It was used by Parker and Raval in their effective action. Therefore, 
it allows us to directly compare our results to the ones that were obtained from their analytic 
solution that neglected derivatives of curvature invariants during the transition to an expan- 
sion of the universe in which R is constant. This choice for 7^2 is sufficiently complicated 
to be representative of other possibilities that contain higher derivatives of the metric. In 
our numerical integration, we do not neglect higher derivatives of the metric in the relevant 
terms, and we carefully analyze the possibility of runaway solutions. As it would be too 
lengthy to giw the details here, we will present the full numerical analysis in a paper now m 
preparation j^]; but here we give our main conclusions and plot a representative numerical 
solution. 

As we noted earlier, there are other possible choices of 7^2- Our numerical analysis using 
the expression in Eq. ()5.9p suggests that the transition from a classical expanding universe 
to an expansion with R constant depends only on certain general features of 7^2, and thus 
may be a fairly robust generic feature of a wide class of asymptotic forms of the heat kernel. 
Therefore, it seems reasonable to also look at effective actions suggested by other underlying 
theories that have similar higher derivative terms. For example, string theory leads to low- 
energy effective actions that contain higher derivative terms analogous to those that appear 
in our low-energy effective action. Thus, our methods and analysis may be of interest to 
string theorists. 

Combining Eqs. ()5.8j) . ()5.5|1 . and ()5.9|1 . and using the variations presented in Appendix 1X1 
[see Eqs. ()A1|) - ()A6|) ]. we have (after dropping the subscripts of X2 and M2 for simplicity): 

(5.10) 

where we have defined 

= VM^ + e2 = v/(m2 + xR)^ + e\ (5.11) 
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and the regular tensors AjJ^, j = 0,l,2,3,4, are given by 

-AR^'^'^R^p - 4/?'^"^^/?"„^^ + m^'^'R^ } , (5.12) 

(1) 360 VVM4 + e2/ ^ a „ a 

+2 (V^V^ - g^'^n - R^"') [2(1 - 30OaR + R^^^^R^fs^s - R'^^Ras] + IR^'^UR 

-60^V^ RVR + 2g'"'R°^VaVpR - 4i?°(^V^Va^ - 8/?^"'''^V«V^i?} (5.13) 

+ [V^VR - gf^'DR + 2V(^/? V'') - 2^'^'^V^i? Va] [2(1 - 30^)ni? + R^^^^^R^/s^s 
-R^f^Rap] + {R^'g'"^ + R'^^g^'' - 2i?°(^^")^ - 4i?'^"'^^) (V„i? V/^i?)} , (5.14) 

1 / /If 6 _ o 2 /\^2 \ 

+ {V^'RVR - g^'^V^RVxR) [2(1 - 30f)ni? + R'^'^^^R^f^^s - R'^'^Rap] 
+2(1 -30f)V"/2V„ {VRVR- g^'^'V^RVxR)} , (5.15) 

= - ^ ( (M^ + e^)^ J " (V^i? V./?) (V-i? V^i? - g^^V-RV^R) . 

(5.16) 

The term (T/^g) in Eq. ()5.1()p stands for the part of the expectation value of the energy- 
momentum tensor coming from the regular part of the effective action. 

Notice that the ratios involving M and e appearing in parenthesis in each one of the 
Eqs. ()5.13|) - ()5.16p are bounded functions of M. In fact, for e very small, these ratios are 
of order 1 except for very particular values of M for which the ratios become close to zero. 
Notice also that in the limit e — all these ratios become equal to 1. In the following 
section we will apply the expectation value given in Eq. ()5.1()j) to an FRW spacetime and 
analyze its cosmological consequences. 



VI. QUANTUM SCALAR FIELD IN AN FRW UNIVERSE: THE VCDM COS- 
MOLOGICAL MODEL 

Now we shall apply the results obtained in the previous sections to a cosmological space- 
time. The goal is to show that the present accelerating expansion of the universe may be 
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explained as due to the nonperturbative-infrared form of the effective action calculated in 
Sec. Such a model for the "dark energy" was first proposed by Parker and Raval and it 
is known as the VCDM or vacuum metamorphosis model jist-jrn]. 

As an idealization, we will consider our universe as being described by a spatially-flat 
FRW spacetime, with line element 

ds"^ = -dt^ + a{tf [dx\ + dxl + dxf) , (6.1) 

filled with non-interacting matter, radiation, and a scalar field with zero expectation value, 
(0) = 0, and a expectation value for its energy-momentum tensor {T^") given by Eqs. ()5.10|) - 
()5.16|) . The semi-classical Einstein equations (with zero cosmological constant) then read 

R"" - \ g^^R = SttGn {T^ + Tl^^ + (T^^)) , (6.2) 

where T^'^ = p.,. {g^'^ + ^u^u'^) /"i and = pmU^u'^, with being the future-pointing, 
normalized vector field orthogonal to the homogeneous and isotropic hypersurfaces, and pr 
and pm being the radiation and matter energy densities, respectively, as measured by the 
family of geodesic observers with four- velocity (comoving observers). Moreover, p^ and 
Pm are constant over each homogeneous and isotropic hypersurface. 

The symmetries of the FRW spacetime greatly simplify the task of solving Eq. (j6.2j) . 
It is not difficult to see that the most general form of a rank-two tensor field T'^'^ that is 
consistent with the FRW symmetries is 

T^"" = Au^u" + 8 {g^"" + u^u") , (6.3) 

with A = T^^^u^Uy and B = (T^ + A) /3 being scalar functions that are constant over each 
homogeneous and isotropic hypersurface, i.e., 

V^^ = -u^V^A and V^i3 = -u^u^V^B . (6.4) 

Take, then, T^*^ to be the difference between the left-hand-side (l.h.s.) and the right-hand- 
side (r.h.s.) of Eq. ()6.2|) . which in an FRW spacetime clearly exhibits the form presented in 
Eq. ()6.3p . This shows that the usually ten (semi-classical) Einstein equations, T^"^ = 0, are 
reduced to two equations, ^ = i3 = 0, in an FRW spacetime. Moreover, considering that 
V^T^'' = is also satisfied [for both sides of Eq. ()6.2|1 satisfy this condition], we have, after 
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using Eq. ()6.4p and the fact that u'^ is a geodesic field, 

= V^T^'^ = u-' [V,, {Au'^) + SV^nl , (6.5) 

which tells us that equation W ^T^^^ = is equivalent to i3 = (VaM°)~^ {Au^), provided 
VaW"" 7^ 0. In this case, ^ = and V^T^'^ = imply that i3 = 0. Then, summarizing 
what we have shown so far, solving Eq. ()6.2|) in an FRW spacetime is equivalent to solving 
A = B = which, in turn, is equivalent to solving ^ = and V ^T^" = (assuming 
VqM" 7^ 0). In other words, we are left with the problem of solving 

R'^'u^u, + Ir = SvrG^ (pr + Pm + {Tn u^u,) , (6.6) 

V^(Tr + Tr) = 0. (6.7) 

[To obtain Eq. (jHHI) we have used the fact that V ^.{R^"' ~ g^"" R/2) = {Tf"") = are iden- 
tities since both yf—g^R^" — g^'^R/2) and y/—g {T^") are given by the functional derivative 
of an action with respect to the metric] Eq. ()6.7p is simple to solve analytically, but the 
Eq. ()6.6p can only be solved numerically. 



Space does not permit us to go into the details [37[ of solving numerically the ordinary 
differential equation for a(t) obtained from Eq. (j6.6p . Here, we show representative plots, 
and summarize the main conclusions of the numerical calculations. In Fig. ^ we show the 
result for the scale parameter a{t) for a universe with present value of matter and radiation 
energy densities such that VLmi) = Pm{to)/Pc = 0.34 and Qro = Pr(to)/Pc = 8.33 x 10~^, 
with pc = 3Hq/ (SttGn), Hq being the present value of the Hubble constant. We compare 
this result (solid line) with the approximation for a{t) given by the earlier version of the 
VCDM model, where a constant-scalar curvature stage follows the usual matter dominated 
phase of the universe (dashed line). We can see that the "constant-i? approximation" in fact 
provides a very good (analytical) approximation to the numerical solution a{t). The mass 
of the scalar field was chosen to be such that rh/Ho = 3.26, and in order to facilitate the 
numerical analysis a non-zero value for e was assumed (see discussion at the end of Sec. IIVI 
for the physical meaning of e). Also for numerical reasons, the value of the dimensionless 
parameter = oGnHq, on which a{t) depends through the quantum energy-momentum 
tensor {T^'^), was taken to be much larger than its physical value ao ~ 10~^^^. This, however, 
should not be a problem since we find that the smaller the value of > 0, the closer the 
numerical solution for a{t) becomes to the constant-i? approximation. 
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Numerical solution (ao=10 £^=W 

Constant-i? approximation 

h ; - 

0.2 0.4 0.6 0.8 1 

Hot 

FIG. 1: Plot of the numerical solution (solid line) and constant- i? approximation (dashed line) for 
the scale parameter a{t), as a function of time t, of a universe with matter and radiation content 
given by i^mo = 0.34 and ri^o = 8.33 x 10^^. In order to facilitate the numerical analysis, we have 
used eno = ^/Hq = 10^^ and uq = uGnHq = 10^^. Also, we used m/HQ = 3.26. 

In Fig. 121 we plot the Hubble parameter H{z) (normalized by its present value Hq) as a 
function of red-shift z. Note again that the numerical solution for H{z) (solid line) is very 
well approximated by the constant-i? approximation (dashed line). 

Finally, in Fig. El we present the numerical evolution of the universe (solid line) in a 
diagram showing its scalar curvature R and the square of the Hubble parameter H^. Notice 
that after spending some time in the matter-dominated stage (represented by the dotted 
line), the universe makes a transition (not as sharp as in the constant-i? approximation) to 
an era dominated by the expectation value of the energy-momentum tensor of the quantized 
scalar field. During this latter era, the universe enters a period of accelerating expansion 
that lasts forever (in the spatially-flat FRW case) and that approaches, asymptotically in 
the future, an exponentially fast expansion with Hubble parameter H m/{2\^). 

It is important to mention that in numerically solving the semi-classical Einstein equations 
presented here, we have taken into account the higher derivative terms. Note that Eq. ()6.6|) 
leads to a fifth-order differential equation for a{t). For this reason, it is natural to expect 
that there exist solutions to Eq. ()6.6|) that are not physically acceptable (e.g., runaway 
solutions). Notwithstanding, we find that taking initial conditions such that a{t) behaves 
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FIG. 2: Plot of the numerical solution (solid line) and constant- i? approximation (dashed line) for 
the Hubble parameter H{z), as a function of redshift z, of a universe with matter and radiation 
content given by i^mo = 0.34 and QrO = 8.33 x 10~^. In order to facilitate the numerical analysis, 
we have used eno = ^/Hq = 10~^ and oq = oGnHq = 10~^. Also, we used fh/Ho = 3.26. 

classically (e.g., describing a matter-dominated universe) during some period of time in the 
past is enough to select only physically acceptable solutions, all of them evolving, eventually, 
to a phase of accelerating expansion. Since the vacuum energy density and pressure are 
negligible at early times, the classical initial conditions are the natural ones to impose. 
Initial conditions sufficiently far away from the classical ones in the past do appear to give 
rise to unphysical solutions, but it is remarkable that the class of natural initial conditions 
evidently gives only physically reasonable solutions. 

As noted earlier, we also numerically integrated the semiclassical Einstein equations that 
result from using the ansatz of Eq. ()3.14j) with n = 0. These equations are simpler to solve 
numerically because they contain fewer time derivatives of the metric. We found that, as 
expected, the rate of change of the scalar curvature R decreases as R approaches the value 

from above. However, instead of making a transition to an expansion in which R remains 
close to m^, the scalar curvature effectively bounces off the value and evolves toward 
increasing values of R. The figures that we show are all for the case of n = 2, which seems 
to be of more physical interest. 
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FIG. 3: Plot of the evolution of the square of the Hubble parameter and scalar curvature R 
of a universe with matter and radiation content given by i^mo = 0.34 and ilro = 8.33 x 10~^. 
Note that after a classical period of expansion (represented by the dotted line), the universe enters 
an era dominated by the energy-momentum tensor of the quantum scalar field that prevents the 
scalar curvature R from dropping below the value m?. Eventually, the universe enters a stage 
of accelerating expansion that leads to an asymptotic exponentially- fast expansion. In order to 
facilitate the numerical analysis, we have used eno = ^/^o ~ 10"^ ctQ = oiGnHq = 10~^. 
Also, we used tti/Hq = 3.26. 

VII. DISCUSSION 

We have reconsidered the theory of the vacuum metamorphosis transition that occurs in 
the vacuum cold dark matter cosmological model from a manifestly nonperturbative point of 
view. We showed first that the terms in the vacuum energy-momentum tensor that become 
large when R is close to the value = —rp?jx derive from the large-s asymptotic behavior 
of the heat kernel. Then by examining the large-s asymptotic form of the exact heat kernel 
in the de Sitter, Einstein static, and linearly expanding FRW universes, we arrived at a 
reasonable ansatz for the dominant asymptotic behavior of the heat kernel in a general 
FRW universe. The key feature of the heat kernel that ultimately causes the vacuum energy 
density and pressure to become large is the presence of a term of the form exp{—ixRs) in the 
large-s regime of the heat kernel. Our approach is manifestly non-perturbative in s because 
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it involves the large-s asymptotic form of the heat kernel, and does not involve summation 
of a series in powers of s. [There is also a factor exp(— im^s) in the heat kernel that is exact, 
and is thus valid for any nonzero mass m and for all values of s.] 

We then obtained the explicit renormalized expression for the terms in the vacuum ex- 
pectation value of Tfj^iy that become large when the scalar curvature R is close to m^. These 
terms involve up to fifth-order time derivatives of the FRW scale factor a{t). We also in- 
troduced a small parameter e (of dimension m^) that phenomenologically describes possible 
weakly dissipative effects, and softens the growth of T^,^ as it! approaches m^. 

Finally, we adopted a specific form for the invariant factor, quadratic in the Riemann 
tensor, that multiplies the exponentials in the heat kernel, and we numerically integrated 
the Einstein equations in a spatially flat FRW universe, with the source consisting of the 
energy-momentum tensor of classical matter and radiation and the expectation value of T^j, 
of the quantized scalar field (f). We numerically evolved the terms having the higher time 
derivatives of a{t). We found that if the universe in the past evolved classically (when the 
vacuum expectation values of T^i, are negligible), then there are no runaway or unphysical 
solutions. All solutions with such classical initial conditions undergo a transition to an 
accelerating expansion that approaches the de Sitter expansion at late times. 

We also noted (see the Introduction) that the theory considered here may provide a 
mechanism for early infiation, which could be caused by heavy bosons, or by the very same 
ultrahght acceletron that may be responsible for the presently observed acceleration of the 
universe. 
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APPENDIX A: VARIATIONS 



We present here the hst of the variations that were used in Sec. |V] (recall that 
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Notice that one can obtain the variation of the factors multiplying the logarithmic term in 
the l.h.s. of the Eqs. ()A1|) - ()A6|1 as the particular case where x = 0. 
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